#include <eigen3/Eigen/Dense>
#include <iostream>
int main(){
    Eigen::MatrixXd matA(3,3);
    matA << 1, 3, 7,
            2, 5, 1,
            1, 8, 9;
    Eigen::MatrixXd matB(3,3);
    matB << 1, 0, 4,
            2, 7, 1,
            1, 2, 3;
    Eigen::MatrixXd matC = matA + matB;
    std::cout << "A + B = \n" << matC << std::endl;
    return 0;
}